Probing protoneutron star density profile from neutrino signals. 
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<N : ABSTRACT 

Oh' Supernovae of Type II is a phenomenon that occurs at the end of evolution of 

^ ! massive stars when the iron core of the star exceeds a mass limit. After collapse 

j^ I of the core under gravity the shock wave alone does not succeed in expelling 

the mass of the star and in this sense the role of neutrinos is the most important 

K^ \ mechanism to do so. During the emission of neutrinos flavor conversion is possible, 

"x" ■ related the phenomenon of oscillations, which however depends directly on the 

^ . particular density proflle of the medium. In this paper we present results of 

O ' numerical simulations of neutrino flavor conversion in protoneutron stars and 

oo ; 

O ■ after collapse. The probabilities of survival for a given flavor in a complete three- 

• j-{ ! flavors framework is estimated through an algorithm which conserves unitarity to 

Vh ' a high degree of accuracy. The sensitivity of the results to the different adopted 

models for the protoneutron star structure is examined in detail demonstrating 

how the neutrino signal could be used to check the validity of models. 

PACS : 21.65.+f , 24.10. Cn , 26.60. +c , 03.75.Ss 



I. INTRODUCTION 



The solution of the solar neutrino "puzzle" [ll, y] has pointed out the relevance 
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of the MSW mechanism [3| in the neutrino propagation through the dense medium 
of astrophysical objects. The mechanism is a consequence of the neutrino masses 
and the corresponding flavor mixing. If neutrino are emitted from the core of 



a star, the initial flavor distribution can be distorted by the flavor conversion 
that can occur during the propagation outwards. On one hand this can be a 
complication in interpreting the possible neutrino signals, on the other hand 
the presence of conversion can be of great value in constraining the models of 
the astrophysical objects. The solar neutrino problem is a typical case where 
flavor conversion offers a solution of the puzzle and at the same time put serious 
constrains on the modeling of the internal structure of the star, which would not 
be easily accessible otherwise. All that has stimulated a large number of studies of 
neutrino propagation in more compact stars, like supernovae, protoneutron stars 



and neutron stars 



fl. 



These astrophysical objects have a much larger central 
density and therefore neutrino flavor conversion can occur only in the outer part, 
where the density can decrease to values compatible with the MSW mechanism. 

Of particular interest is the case of protoneutron stars, where the neutrinos of 
all three flavors can be considered in thermal equilibrium and thermally generated 
at the center of the stellar object. In principle the neutrino energy spectrum 
of each flavor is determined by the temperature at which they decouple from 
matter, i.e. the extension of the corresponding neutrino sphere. However, due to 
possible conversion through the MSW mechanism, the energy spectrum and each 
flavor fraction is modified after the propagation through the stellar medium. This 
modification depends crucially on the density profile, and therefore the neutrino 
signal can be used in principle to constraint the models of the protoneutron star 
structure and evolution. The effect of the shock wave on the neutrino signal 
has bee. analysed ,. seve... papers fl fl Q. both as a fu.et.o. of t„„e and ,„ 
relation to neutrino mixing parameters. The dependence on the structure of the 
progenitor star has been analyzed in ref. [8|, while the expected relic neutrinos 
characteristics have been studied in ref. [9|]. 

However, to the authors knowledge, no systematic study of the relation be- 
tween possible neutrino signals and protoneutron star structure have been carried 
out in the literature. In this work we present a detailed analysis of the neutrino 
propagation in the protoneutron star medium and we try to relate directly some 
characteristics of the neutrino signals to the star density profile, which, in prin- 
ciple, could be used to distinguish among different models or model parameters. 

In this study we will neglect neutrino-neutrino interaction, which has been 



recently shown [5] to be able to produce giant coherent flavor conversion. This 
coherent process can occur under particular conditions, and, to simplify the anal- 
ysis, we do not consider this possibility in the paper. 

II. GENERALITIES AND NUMERICAL METHOD. 

It is customary to schematize the propagation of neutrinos outside the neu- 
trino sphere as coherent and subject only to a single particle potential, i.e. only 
forward scattering is present [lO(] . The medium has then only the effect of shift- 
ing the neutrino masses locally at each point of the star. As it is well known, 
due to the dominance of electrons in the lepton content of the protoneutron star, 
this shift is flavor dependent and the neutrino mass spectrum is strongly altered. 
In fact, due to the electron dominance, electron neutrinos interact also through 
charged weak current, while the r and /i neutrinos interact essentially through 
neutral weak current only. In the deep interior of the star the electron neutrino 
mass is shifted at much higher value then the neutrinos of the other two flavors. 
The mass spectrum is then formed by a singlet eigenstate, corresponding to pure 
electron neutrinos, and an eigenstate doublet below, involving the other two fla- 
vors. During the neutrino propagation outwards, the electron density decrease 
and the mass difference between the singlet and the doublet decreases, but to the 
extent that the mass difference stays large enough, the structure of the spectrum 
remains unaltered. Since, at least in the so called "direct hierarchy" scheme, the 
lowest neutrino mass in free space has a large electron flavor content, at a certain 
electron density the singlet eigenstate must cross the doublet and strong mix- 
ing can occur, i.e. neutrinos can jump among the different local (instantaneous) 
mass eigenstates, thus altering the neutrino populations among them. This is the 
essence of the MSW mechanism. According to the degree of adiabaticity of the 
process, some flavor conversion can occur, which determines the flnal flavor dis- 
tribution. If the crossing is fast enough, the neutrinos in the initial electron flavor 
singlet eigenstate will jump in the lowest eigenstate, which will flnally merges in 
the lowest neutrino mass in free space. The latter has a deflnite flavor content 
(mainly a mixture of electron and muon flavors). On the contrary, if the process 
is completely adiabatic, the neutrinos will remain in the same (higher) local mass 



eigenstate, and finally they will propagate in the highest mass eigenstate in free 
space (mainly a mixture of r and fi flavors). 

From this simple considerations it is clear that the crucial parameter of the 
process is the adibaticity and therefore the density profile, which fixes the local 
density slope, has a direct influence on the final neutrino signal. It is this connec- 
tion between the density profile and the final mass and flavor populations that 
we want to exploit systematically in this paper. 

In a simplified description where only two mass eigenstates are considered most 
of the calculations cannot be done analytically, and a fortiori in a more realistic 
description with all three flavors included only numerical simulations are possible. 
Following the standard treatment, the initial flavor distribution can be evolved 
by a 3 X 3 matrix. In the case of three neutrino the flavor eigenstate Ue, v^j, and 
Vj. and mass eigenstate vi, v^ and 1/3 are related through the linear relation 
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In general, in the case of Dirac neutrinos the lepton mixing matrix U in eq. ([T]) 
depends on three mixing angles ^12, ^13 and 6'23 and one CP- violating phase 5, 
while in the case of Majorana neutrinos there are two additional, so-called Ma- 
jorana phases. It is convenient to use the parametrization of the matrix U which 
coincides with the standard parametrization of the quark mixing matrix 
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sin^jj. The probabilities of oscillations between 



various flavor states, unlike in the two-flavor case, are in general do not have a 
simple form. The evolution equation describing neutrino oscillations in matter is 
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Here U is the 3x3 flavor vacuum mixing matrix defined in Eq. ([T]) and 
Vcc{r) is the matter- induced neutrino potential for unpolarized medium of zero 
total momentum defined in terms of electron number density Ne{r) and Fermi 
constant Gp as 

V^r) = V2GFN,{r). (4) 

It is rather difficult to study this equation analytically, even in an adiabatic 
approximation. As we will see, the numerical treatment is necessary also because 
in some cases the density profile does not produce a real MSW resonance, but 
rather a "modulation" of the flavor populations, which can be hardly treated an- 
alytically. This requires a quite stable numerical algorithm in solving iteratively 
eq. ([3]). We used a method, based on the Hamilton-Caley theorem, which guaran- 
tees the fulfillment of unitary to a very high degree of accuracy and it is therefore 
suitable for calculations which extend for large time an long distances. Typically 
the calculations were performed from the center to the stars up to 10^ kilometers. 

Contrary to other methods used to solve the evolution equations of the neutrino 
states by means of techniques like Runge-Kutta or similar our numerical algorithm 
using the Hamilton-Caley's theorem assures unitarity of the evolution operator 
to a very high degree of accuracy. If the wave function of the flavor triplet of 
a neutrino is written as a three components vector tpi^ = (z/g, z/^, z/^) the flavor 
evolutions is completely determined by the Schrodinger equation in the compact 
form 
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where the matrix H is the Hamiltonian operator. The code obtain ip^ll + 61) 
from ipu{l) using an iterative procedure step by step: 
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preserving the unitary of ipi, automatically and becoming exact if H is indipen- 
dent of spatial coordinate. Therefore the step sizes employed in numerical codes 
are not restricted by the size of H but are restricted by the rate of change of H. 
In this algorithm the Hamilton-Caley's theorem is used to calculate step by step 
the exponential operator in the eq. i^ to determinate the evolution of ip^, both 
in the simple two flavors framework and in the standard three flavor framework. 

III. NEUTRINO SIGNALS AND THE PROTONEUTRON STAR 
STRUCTURE. 

Numerical simulations of the collapse and evolution of supernovae up to the 
stage of protoneutron star is quite difficult and affected by many uncertainties. 
For these reason simulations up to few tens of seconds after collapse are quite rare. 
The reference simulations are those of the Livermore group ll| , which are often 
used in the literature for different types of analysis. These simulations include 
also the shock wave front. Phenomenological models are also used for qualitative 
considerations [a]. 

At the protoneutron star stage of the evolution neutrinos are essentially ther- 
mally produced. If all three flavors had the same temperature at the moment of 
decoupling from matter, i.e. the neutrino sphere is independent of flavor, then all 
three flavors had the same fraction of 1/3 and the same energy spectrum. Any 
possible conversion mechanism would not affect such a population distribution 
and no characteristic neutrino signal could be present. It is therefore essential 
to specify the initial temperature and chemical potential for each flavor in the 



core of the star. For deflniteness we will follow ref. 12j, but with the possibility 
of slightly modify the parameters to check the sensitivity of the results on their 
speciflc values. Each flavor population is described by a Fermi distribution 

where rji, is the degeneracy parameter, T^, is the neutrino temperature and 

x'^dx 
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is the Fermi integral with rank k = 2. 

The values of the degeneracy parameters are taken r]^^ = rj^^ = rj^^ = ?> and 
neutrino temperatures are taken T^^ ^ 2.76 MeV and T^ = T^^ ^ 6.26 MeV. 
The normalization factors F2{ri^) specify also the total content of each flavor. 
Once the initial conditions are fixed, one has to define the parameters of the 
flavor mixing 3x3 matrix. Presently the only parameter that remain uncertain 
is the so-called 613 mixing angle, while for all the others, standard established 
values can be used [J]. We choose the value sin^ 6*13 = 10~^. 

Since the density profile changes with time also the final neutrino signal is 
expected to change with time. For illustration, at the fixed time of t = 4 s 
and for the energy E = 3 MeV the flavor survival probabilities are reported 
in Fig. [1] as a function of the distance from the core of the star up to radius of 10^ 
kilometers, where the influence of matter is surely not relevant. For this reason 
the probabilities or transition probabilities are usually rapidly oscillating in the 
final part of the propagation. We have then followed an averaging procedure over 
many oscillation period to establish the final average value of each flavor content. 
This is justified by the wide region where neutrino are actually emitted, which 
introduces a random phase. In any case, one can easily recognize the place where 
the MSW resonances occur. 

This actually corresponds to the higher resonance, i.e. the place where the 
initial highest mass eigenvalue crosses the first one below. 

In each row is reported the survival probability of each flavor (along the diago- 
nal) and the contribution to each flavor coming from the other two flavors (off di- 
agonal panels) . In other words each panels contains the evolution of the flavor con- 
version probabilities along the neutrino path. At each radial distance r the sum of 
the values reported in a given row gives the total survival probability of each fla- 
vor. It has to be noticed the relevance of the "re-population" of e.g. the electron 
flavor by the other two flavors, without which the electron flavor content would 
drop to zero. Another example with different values of energy and fixed time is in 
Fig. [21 Of course the density profile changes during the neutrino propagation but 
the main conversion processes occur in a time interval short enough that the den- 
sity profile can be considered fixed. In this sense the nominal value of time t{s) 
must be taken just as a nominal value labeling the particular case considered. 



From figures 1 and 2 it is apparent the interplay of tlie MSW conversion 
process and tlie non-resonant modulation due to tlie density profile. This is 
mainly due to the fact that the MSW resonance width is quite wide, and some 
degree of flavor conversion can occur even if the resonance density is not strictly 
crossed. The final signal depends in a complex way on the density profile. This 
also indicates the need of numerical calculations. It has to be stressed that the 
density at which the MSW resonance is expected to occur is quite low with respect 
to the typical nuclear density and in this density region matter is dominated by 
electrons, positrons and photons. However, the overall density profile depends 
also indirectly on the nuclear Equation of State (EOS), since the EOS is one of 
the main physical ingredient which determines the evolution and structure of the 
protoneutron star, and thus the density profile even in its outer part. 

We have repeated these calculations at different times and different neutrino 
energies. The final neutrino energy spectrum of each flavor can be then extracted 
at different protoneutron star evolution time. This can be obtained by multiplying 
the initial Fermi distributions by the survival probability. In Fig. [3] is reported 
the fraction of electron neutrino for two different times as a function of neutrino 
energy. It is clear that in these cases the electron neutrino fraction is only slightly 
modulated as a function of time. Also the initial spectrum is affected by the 
propagation through the protoneutron star matter, however the final spectrum is 
expected to change only slightly at different times. 

Indeed the spectrum is reported in Fig. (jlj) at the two different times in 
comparison with the initial one. It has to be noticed the strong distortion of 
the initial spectrum, due to the conversion process, but this distortion depends 
little on time. Of course the overall normalization depends on the total neutrino 
flux, but the considered properties, i.e. the electron neutrino fraction and the 
corresponding spectrum, do not depend on the value of the flux. 

To see the sensitivity of these results on the density profile we have repeated 
the calculations by changing artificially the density profile of the protoneutron 
stars. To this purpose we have followed the profile parametrization introduced 
in ref. |6j, according to which the density profile can be effectively approximated 
for post-bounce times t < 1 s as 
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and for post-bounce time t > 1 s as 



f ■ fir) , r ^tq , 
p{r) = po(r) • r (10) 

1 , r^rs , 

where ^ ~ 10 and the function /(r) parametrizes the rarefaction zone above 
the shock front. The parametrization for /(r) is taken 

In /(r) = [0.28 - 0.69 In(r5/Km)] ■ [arcsin(l - r/rg)]^-^. (11) 

In eq. ( TTOl) and eq. ( ITTl) it is assumed a shghtly accelerating shock-front 
position xs, 



rs{t) = r° + vst + -ast" (12) 

with parameters fixed by 



Xs — —4.6 X lO^km , 

t;s~ 11.3 X lO^km/s , (13) 

as ~ 0.2 X 10^km/s^ 
The value of the parameter a = 2.4 corresponds to the parametrization which 



reproduce closely the simulations of ref. [ll|, and this is the value used in the 
calculations reported above. We then varied the value of a to obtain a steeper 
or a more shallow density profile, which should possibly reflect different model 
simulations of the proton-neutron star dynamics. This is illustrated in Fig. [5] , 
where the density profiles for different values of a are compared at fixed times. 

The modification of the profile has a direct effect on the neutrino fraction 
as shown in Fig. [6] for the case of a = 2.1. At the lower energy, below about 
15 MeV, the electron neutrino fraction presents a sharp drop at an early time 
evolution of the protoneutron star. This means that the electron neutrino energy 
spectrum should be drastically suppressed in the low energy part. At a later stage 
the energy spectrum should change again. The final spectra at the two different 
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times, in comparison with the initial Fermi-Dirac spectrum, for electron neutrinos 
are reported in Fig. [71 Also in this case the final spectrum is strongly distorted 
by the conversion process which occurs during the neutrino propagation in the 
protoneutron star matter, but what is relevant is that the distortion depends 
now more clearly on time. The increase of the high energy tail is due to the 
conversion of fi and r neutrinos in electron neutrinos. Modification of the initial 
Fermi-Dirac distribution can also occur due to energy dependence of the neutrino 
transparency [l3| , but the type of dynamical evolution of the spectrum reported 
in Fig. 7 should be a distinct feature of a shallow density profile, also independent 
of the total neutrino fiux, which also depends on time. Comparison between Figs. 
IS1B1 and Figs. [^|71 shows the sensitivity on the density profile of the protoneutron 
stars during its evolution. 

The reason for such a behavior is directly linked to the higher degree of adi- 
abaticity of the MSW mechanism, in the wider sense discussed above. In fact, 
increasing the value of a to 2.7 both electron neutrino fraction and spectrum 
become almost independent of time, see Figs. 8,9. 

We have also tested the dependence of the results on the assumed values of the 
neutrino temperatures. Decreasing the temperature differences all the considered 
effects decrease, which indicate the possibility of measuring the original neutrino 
temperatures. 

IV. CONCLUSIONS. 

In this paper we have presented evidence of possible neutrino signals from pro- 
toneutron stars based on the flavor conversion processes which can occur inside 
the star and which are sensitive to the density proflle of these astrophysical com- 
pact objects. Generally speaking a more shallow density proflle would produce a 
percentage of electron neutrino strongly suppressed at low energy, approximately 
below 15 MeV, few second after collapse. This suppression should disappear after 
10-20 seconds. Correspondingly, the electron neutrino spectrum should show a 
characteristic time dependence. 

Of course, all these signals require an appreciable number of events in the 
detectors, and only with the next generation of neutrinos detectors one can hope 
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to have indications of the signals studied in the present paper. 
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FIG. 1: Survival probabilities and transition probabilities for different flavors along 
neutrino path inside the protoneutron star matter. The energy of neutrinos is 3 MeV 
and the time after bounce is t = 4 s. The density profile is characterized by the 
parameter a = 2.1, see the text for detail 
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FIG. 2: The same as in Fig. 1 but at the different time t = 12 s. 
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FIG. 3: Fraction of electron neutrinos as a function of energy at t = 4 s (full line) and 
t = 27 s (dashed line) after bounce. The parameter of the profiles has been fixed at a 
= 2.4 . 
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FIG. 4: electron neutrino spectrum corresponding to the case of Fig. 3. The full line 
is the initial Fermi-Dirac spectrum, the dashed line is the spectrum at t = 4 s and the 
dotted line at t = 27 s. 
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FIG. 5: Density profiles at different times after bounce corresponding to a = 2.1 (left) 
and to a = 2.7 (right). 
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FIG. 6: The same as in Fig. 3 but for a = 2.1 
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FIG. 7: The same as in Fig. 4 but foe a = 2.1 . 
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FIG. 8: The same as in Fig. 3 but for a = 2.7 
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FIG. 9: The same as in Fig. 4 but for a = 2.7 . 
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